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Abstract: Autophagy plays an important role in cancer and it has been suggested that it functions not only as a tumor 
suppressor pathway to prevent tumor initiation, but also as a pro-survival pathway that helps tumor cells endure metabolic 
stress and resist death triggered by chemo therapeutic agents, including acquired resistance. We aimed to identify small- 
molecule autophagy inhibitors using a HTS/HCA approach through a phenotypic, cell image-based assay, in order to 
screen multiple biological targets simultaneously and to screen compounds in a physiologically relevant environment. 
LC3 is a component of the autophagosome, which undergoes a cytoplasmic redistribution from diffuse to punctate dots 
during autophagy. We employed HeLa cells stably expressing EGFP-LC3 in a primary phenotypic screen. As a first step, 
a "Validation Library" of about 8,000 pre-selected compounds, about 25% of which had known biological activity and the 
others representing a range of chemical structures, was run in duplicate both to assess screening suitability and likely hit 
rate, and to give a valuable preview of possible active structures or biological targets. The primary screen of about 0.25 
million compounds yielded around 10,500 positive compounds. These were tested in a suite of further cellular assays de- 
signed to eliminate unwanted positives, together with the apphcation of chemi- and bioinformatics to pick out compounds 
with known biological activity. These processes enabled the selection of compounds that were the most promisingly ac- 
tive and specific. The screening "tree" identified, amongst others with as yet unidentified targets, chemical series active 
against autophagy-relevant biological targets ULK or Vsp34, validating the phenotypic screening methods selected. Fi- 
nally, about 400 compounds were fully qualified after following this triage. The development of the assays, compound 
screening process and the compound triage is described. 
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INTRODUCTION 

Autophagy (from the Greek, "auto": oneself; 
"phagy": to eat) refers to a cellular degradation pathway that 
involves delivery of cytoplasmic cargo sequestered inside 
double-membraned vesicles to the lysosome. The cellular 
events during autophagy follow distinct stages (Fig. 1): vesi- 
cle nucleation (formation of the isolation mem- 
brane/phagophore), vesicle elongation and completion, fu- 
sion of the double-membraned auto phagosome with the 
lysosome to form an auto-lysosome, and finally, lysis of the 
autolysosomal contents and recycling of nutrients. Basal 
autophagy is necessary for most cells to eliminate protein 
aggregates and damaged organelles and can be rapidly 
upregulated during various pathological and physiological 
states, such as starvation. Starvation-induced autophagy en- 
ables the degradation of proteins, carbohydrates and lipids, 
which allows the cell to adapt its metabolism and meet its 



*Address correspondence to these authors at the Sanofi Tucson Innovation 

Center, 2090 Innovation Park Drive, Tucson, Arizona 85755, U.S.A; 

Tel +1 520 544 6875; Fax +1 520 544 6804; E-mail: jane.peppard@sanofi, 

and Oncology R&D, Sanofi, 94403 Vitry-sur-Seine, France; 

Tel +33 1 58 93 89 93; Fax: +33 1 58 93 34 71; 

E-mail: benoit.pasquier@sanofi.com 



energy needs. Cancer is the first disease that was linked to 
autophagy. The prevailing current view is that autophagy 
functions both as a tumor suppressor pathway that prevents 
tumor initiation, and as a pro-survival pathway that helps 
tumor cells endure metabolic stress and resist death triggered 
by Chemother apeutic agents. Autophagy may also contribute 
to acquired resistance of tumor cells to chemotherapeutic 
agents [1,2]. 

The cellular pathway to autophagy has been well studied 
on the molecular level [2, 3]. Autophagy is negatively regu- 
lated by the serine/threonine kinase mammalian Target of 
Rapamycin (mTOR) which acts upstream of the UNC-51- 
like kinase 1 (ULKl) complex. Autophagy can be induced 
by inhibiting mTOR complex 1 (mTORCl) using starvation 
or inhibitors such as rapamycin. Another key regulator is the 
Beclin 1 complex, a regulatory platform, consisting of Beclin 
1, the Class III phosphatidylinositol 3-kinase (PI3KC3A/^- 
ps34) and the autophagy-related gene (Atg) protein Atgl4L. 
The activity of PI3KC3/Vps34 is regulated during starvation 
and generates the phosphatidylinositol 3 -phosphate 
(PtdIns3P) required for the onset of autophagosomal 
membrane nucleation. Two ubiquitin-like protein conjuga- 
tion systems are also necessary for the formation of the 
autophagosome, the Atg5-Atgl2 and the microtubule- 
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Fig. (1). Schematic outline of the autophagy process. LC3 distribution in the cytoplasm changes once autophagy is initiated as the molecules 
cluster in double membraned organelles called autophagosomes. These fuse with the single membraned lysosomes, and finally the auto- 
phagosomes are degraded. Hydroxychloroquine inhibits the lysosomal degradation step. 



associated protein 1 light-chain 3 (MAP1LC3/Atg8, hereaf- 
ter called LC3) systems. In the LC3 system, cytosolic LC3 is 
cleaved by the protease Atg4 (LC3-I) then converted to 
phosphatidylethanolamine (PE) at the autophagosomal 
membrane (LC3-II) by ATG7/ATG3. LC3 has multiple 
functions in autophagy, including expansion and fusion of 
autophagosomes. LC3 status can serve as a biomarker for 
autophagy, and by fluorescently tagging LC3, auto- 
phagosomes can be detected and their localization patterns 
observed [4, 5]. 

Because of the lack of specific autophagy inhibitors and 
the potential therapeutic applications in cancer [1, 2], we 
aimed to identify small molecule autophagy inhibitors using 
the Sanofi chemical library. We chose to approach the screen 
with a phenotypic, cell image-based assay, rather than select- 
ing from the multitude of potential molecular targets, in or- 
der to cast our net wide and to screen compounds in the more 
physiologically relevant environment of intact cells. In the 
primary screen, HeLa cells stably transfected with EGFP- 
LC3 were starved and the formation of EGFP-containing 
"granules" was quantitated by cell image analysis to identify 
inhibitors of autophagosome formation, an early step in the 
autophagy process. We also designed a screening tree with 
the aim to first identify all possible autophagy inhibitors with 
the primary screen and then actively filter this compound set 
in order to home in on the compounds most specific and 
suitable for our targeted pathway. A suite of phenotypic cel- 
lular assays was set up and only compounds qualifying 
through this filtering process were accepted as true positives. 
Here, we describe the screening progression and the immedi- 
ate outcome of this phenotypic approach to small molecule 
screening. 



MATERIALS AND METHODS 

Reagents & Equipment 

Unless otherwise specified, all chemicals were purchased 
from Sigma-Aldrich (St. Louis, MO); tissue culture reagents 
and supplies from Life Technologies (Carlsbad, CA), Corn- 
ing (Tewksbury, MA) or Thermo-Fisher Scientific 
(Waltham, MA). Draq5 was purchased from Biostatus 
(Shepshed, UK). Multiwell plates were purchased from 
Greiner BioOne (Monroe, CA), Corning or Matrical Bio- 
science (Spokane, WA). Bulk reagents (cells, medium, etc.) 
were dispensed to all plate wells with a Multidrop Combi 
(Thermo-Fisher) or a Matrix Wellmate fitted with plate 
stackers (Thermo-Fisher), and compounds transferred with a 
CyBi®-Well 384-channel simultaneous pipettor (CyBio, Bos- 
ton, MA). Serial dilutions of compounds were performed 
using a BioMek FX® (Beckman Coulter, Brea, CA). Multi- 
well plates were washed with either a Skatron EMBLA-384 
(Molecular Devices, Sunnyvale, CA) or Bio-TekEL405™ 
(Winooski, VT) platewasher; washers were presterilized with 
bleach, then ethanol and flushed thoroughly with sterile wa- 
ter before use. The IN Cell Analyzer 3000, a confocal line- 
scanning instrument, along with cell analysis software mod- 
ules, was purchased from GE Life Sciences (Piscataway, 
NJ). 

HeLa EGFP-LC3 Clone €10 Cells 

EGFP-wild type LC3B plasmid was constructed as fol- 
lows. The LC3B coding sequence was amplified from EST 
clone BC015810.1 (IMAGE) by PCR and then cloned into 
pEGFP-Cl (Clontech, Mountain View, CA). The resulting 



Identification of Autophagy Inhibitors using High-Content Screening 



Current Chemical Genomics and Translational Medicine, 2014, Volume 8 5 



EGFP-LC3 He\s All 




Clone C10 



Cell sorting f 



A 



1} 




EGFP 



EGFP 



1*" 



Fig. (2). Flow cytometry analysis of the EGFP-LC3 HeLa cell population before clonal selection (left panel, blue line) and after (right panel, 
green line), where 8 % of the population were expressing EGFP-LC3, compared to WT HeLa cells (red line). Clone CIO (green line/green 
arrow) was produced by single-cell sorting of the mixed population, expanded, and selected because of its high fluorescent signal. 



plasmid allowed the expression of EGFP-LC3B fusion pro- 
tein under CMV promoter control. HeLa cells (ATCC, Ma- 
nassas, VA) were stably transfected with EGFP- wild type 
LC3B plasmid using Fugene 6 (Roche Applied Science, In- 
dianapolis, IN). Cells were then sorted based on EGFP inten- 
sity using a FACSDiva (Becton Dickinson, San Jose, CA) 
and single-cell derived clonal populations with high EGFP 
expression generated. HeLa EGFP-LC3 Clone CIO was se- 
lected based on fluorescence intensity (Fig. 2) using a FAC- 
SCalibur (Becton Dickinson) and on EGFP-LC3 grains ob- 
served by iCyte (Compucyte, Cambridge, MA). A lack of 
Pgp expression, desirable when screening compounds that 
may otherwise be pumped out of the cell, was confirmed 
using PE-conjugated anti-Pgp antibody (Becton Dickinson) 
with cells analyzed on the FACSCalibur. EGFP expression 
in these cells was monitored over time and found to be stable 
until at least passage 18. In preparation for screening, a large 
single batch of cells at passage 5 was frozen in vials in liquid 
nitrogen, and to maintain assay performance, cells were kept 
in culture only until passage 18 and then replaced with a 
fresh vial of cells. 

Cell Assays Used in the Screening Tree 

Three assays were set up for identifying active inhibitors 
of autophagy. For simplicity, the methods detailed in this 
section are the fmal protocols used in screening. Details of 
the assay development strategy are described in Results. 

1. Primary Screen: Cell Imaging Assay for Autophagy 

HeLa EGFP-LC3 Clone CIO cells were expanded in 
growth medium (DMEM high glucose with glutamine, sup- 
plemented with 10% FBS (heat inactivated, qualified), 100 
units/ml Penicillin/Streptomycin and 1 mg/ml Geneticin) and 
detached from culture flasks using 0.05% Trypsin-EDTA. 
Cells were plated (4,000/well in 50 |il growth medium) into 
384 well black-walled, glass bottom assay plates (Matrical) 
and incubated at 37''C/ 5% CO2 overnight. After washing the 
cells three times (80 |il Earle's Balanced Salt Solution, 
EBSS) with an EMBLA-384 washer, 40 |il EBSS was left in 
the wells and 5 |il compounds or controls added directly at 
lOx fmal concentration. Screening compounds were tested at 
10 |iM/0.3%) DMSO (fmal); autophagy inhibitor control was 
33 nM wortmannin/0.3%) DMSO, negative (starvation, posi- 
tive autophagy) control was 0.3%) DMSO. To inhibit 



lysosomal breakdown of phagosomes, hydroxychloroquine 
sulfate (HC) was diluted in EBSS and 5 |il was added to all 
wells (10 |iM final). Cells were incubated at 37"C/ 5% CO2 
for 2 hr. 

To fix and stain the cells, 25 |il of 12%) formaldehyde in 
Dulbecco's phosphate-buffered saline (DPBS) was added to 
all wells (4%) final). Plates were incubated at 4°C overnight, 
fixative was aspirated from wells and cells washed three 
times with 80 |il EBSS using a plate washer with a final aspi- 
ration leaving 15 |il volume. The nuclear stain Draq5 (50 jil 
in DPBS, 146 nM final) was added and incubated in the dark 
at room temperature for 2 hr whereupon cells were washed 
three times with 80 |il DPBS, leaving the final 80 |il in the 
well. Plates were sealed with black tape and stored at 4°C. 
Plates were imaged within 24 hr on the IN Cell Analyzer 
3000 using 488 nm (excitation) and emission filter 535-45 
nm for EGFP and 647 nm/695-55 nm for Draq5, and ana- 
lyzed in-line using the Granularity Analysis Module GRNl 
image analysis module of the IN Cell analysis software. 

2. Counters creen: Effect of Compound on Cell Viability 
using the Parental HeLa Cell Line in Full Growth Medium 

Cell viability of wild type (WT) HeLa cells was deter- 
mined in the presence of test compounds at the same concen- 
tration used in the primary screen. Cells were detached from 
culture flasks using 0.05%) Trypsin-EDTA and plated into 
Falcon 384-well white solid bottom assay plates: 1,500 cells 
in 45 jil growth medium were added to each well and plates 
incubated at 37''C/ 5% CO2 overnight. Cells were treated 
with 5 |il test compounds at 10 |iM/ 0.3%) DMSO, where 
maximum toxicity (high control) was defined by 3.3 
|iM staurosporine/0.3%) DMSO and compared to the low 
control defined by 0.3%) DMSO. After a further 24 hr incu- 
bation, viability was determined with the CellTiter-Blue® 
Cell Viability Assay (Promega, Madison, WI) as recom- 
mended by the manufacturers, and the fluorescent signal was 
read on an EnVision® plate reader (Perkin Elmer, Waltham, 
MA). 

3. Secondary assay: effect of compound on cell viability in 
an alternate tumor cell line H1299 with or without nutrient 
starvation for 48 hr 

The primary screen was aimed at identifying compounds 
that inhibit autophagy in starved HeLa cells. Ideally for our 
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Fig. (3). Plate map of controls and compounds for primary screening. Compounds were run at 10 jiM in singlicate (320/plate), along with 16 
of each control type on every plate. The results from the compounds were normalized as % plate controls, where appropriate for the feature 
measured; features such as Npass (nucleus number) were treated as raw data. 



purpose, these compounds should reduce cell viability in low 
nutrient conditions but have little or no effect in high nutrient 
conditions. A 24 hr full-growth medium assay of HeLa cells 
was used for the counterscreen but as a final check we added 
a secondary assay into the screening tree to compare in par- 
allel the effect of compounds on cell viability in the presence 
or absence of nutrients. As an added filter a different cell line 
was selected, the human lung carcinoma cell line HI 299. 

HI 299 cells (obtained from ATCC) were expanded in 
growth medium (RPMI-1640 supplemented with 10% FBS 
(Heat inactivated, qualified) and 100 units/ml Penicil- 
lin/Streptomycin) and detached from culture flasks using 
0.05% Trypsin-EDTA. Either 500 (for growth medium) or 
2,000 (for starvation medium) cells in 45 |il growth medium 
were added to each well of 384-well white solid bottom tis- 
sue culture plates using a WellMate and plates were incu- 
bated at ?>TCI 5%) CO2 over night. The next day, medium 
was removed and replaced as appropriate with either 45 |il of 
growth medium or starvation medium (EBSS plus 100 
units/ml Penicillin/Streptomycin). A 5 |il volume of com- 
pound (lOx fmal concentration) in EBSS was added and the 
plates were incubated for a further 48 hr. Cell viability was 
assessed with CellTiter-Blue® Cell Viability Assay as before 
and the fluorescent signal was read on an EnVision® plate 
reader. 

Data Analysis 

Data were analyzed using Microsoft Excel (Redmond, 
WA), Spotfire® (TIBCO, Palo Alto, CA), and GeneData 
Screener® (Lexington, MA), as well as automated image 
analysis using the software provided with the IN Cell as al- 
ready described. Images and associated data were stored on 



the Cellomics'' Store database (Thermo-Fisher) and retrieved 
and further processed for analysis via Screener®. Compound 
plate maps were retrieved from ActivityBase (IDBS, Bridge- 
water, NJ) and associated with image-based data via 
Screener®. In the primary assay, compound plates were set 
up as shown in Fig. (3) and run at 10 |iM in singHcate. Aver- 
aged intra-plate assay controls (maximal inhibition vs. no 
inhibition) were used to normalize results obtained with the 
compounds in each plate, and determine the plate Z' [6]. 

RESULTS 

Cell Assay and Image Analysis to Detect Early Auto- 
phagy: Assay Development 

In the initial set-up for image analysis, HeLa-EGFP-LC3 
Clone 1 0 cells were nutrient starved and then incubated for 2 
hours with test compound in the presence of the lysosomo- 
tropic agent Hydroxychloroquine (HC). Inhibiting lysosome 
function with HC blocks the final degradation of EGFP-LC3 
in the lysosome [4, 5] and thus increases the numbers of 
fluorescently labeled autophago somes available for imaging. 
Cells were imaged on the IN Cell Analyzer 3000, an auto- 
mated, line-scanning confocal cell imager, with on-the-fly 
cell feature analysis using the Granularity Analysis Module 
GRNl. In this algorithm, cells are first identified by their 
nuclear stain and a measurement region for each cell is de- 
fined from the nucleus. Intense fluorescently areas ('grains') 
within this region are identified and analyzed (Fig. 4). In our 
assay, these grains represent the accumulated EGFP-LC3 
labeled autophago somes. The algorithm reports features for 
each cell selected for analysis in the well and well-level sta- 
tistics are computed from these cell-level features. The fea- 
tures selected for analysis are shown in Table 1. 
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Fig. (4). Images of HeLa-EGFP-LC3 cells after 2 hr starvation; cells were stained with DraqS to indicate nuclei (red, Channel 1), EGFP-LC3 
is shown in green (Channel 2). The algorithm Granularity Analysis Module GRNl was apphed to first identify cell cytoplasm by expanding 
out a set number of pixels from the nucleus boundary designated in Channel 1 (large box) and then identify EGFP-autphagosomes located 
within the cell mask, which are counted as "grains" (small boxes). Other measurements are also made, e.g. brightness of objects in each 
Channel; nuclear size and cytoplasm area. Typical images of the assay controls are shown; the left panel shows cells after 2 hr starvation and 
the right shows cells after 2 hr starvation in the presence of the inhibitor wortmannin. 

Table 1. Measurements Selected for Quantitation from the IN Cell GRNl Module 



Feature 


Description 


Ngrains 


Average number grains/cell 


%Pos 


% cells in a well > set threshold [e.g % cells > 6 grains] 


Agrains 


Average fractional grain area/cell [(Total Grain Area/Total Cell Area) x 1000] 


Fgrains 


Average fractional grain fluorescence intensity/cell [(Total Grain Intensity/Total Cell Intensity) x 1000] 


Npass 


Number of cells analyzed 


Inuc 


Average GFP intensity in nuclear region 


Icyt 


Average GFP intensity in cytoplasmic region 


Dmrk 


Average nuclear diameter 



As a first step, optimal concentrations of HC and wort- 
mannin were determined. HeLa-EGFP-LC3 cells were set up 
in 384-well plates as described in IMaterials and IVlethods, 
transferred to BBSS and treated for 2 hr with a concentration 
range of both compounds as indicated. The images were ana- 
lyzed to determine the % cells/well with, initially, >4 grains 
per cell. Incubation with either 5 or 10 [iM HC produced 
similar ICso values for wortmannin (Fig. 5). We chose 33 
nM wortmannin to provide the positive (max. inhibition) 
control. At this concentration of wortmannin, increasing HC 
was found, as expected, to improve the signal and the Z' 
values [6] by inhibiting the dynamic lysosomal degradation 
of the autophago somes: 10 [iM HC was chosen for inclusion 
in the assay. 

Next, the effects on the assay signal of increasing com- 
pound incubation times were compared, adding 10 [iM HC 
immediately after compound addition. The response to 



wortmannin was already maximal by 2 hr starvation (Fig 
6A). Therefore, we selected a starvation time of 2 hr post 
compound and HC addition. 

The IN Cell 3000 is a semi-confocal instrument and until 
this stage we had been using plastic imaging plates, but now 
switched to glass plates from IMatrical, glass being a more 
suitable medium for high resolution imaging. The number of 
grains resolved per cell was observed to increase upon this 
plate change due to the sharper images obtained. We there- 
fore increased the response threshold, from >4 to >6 
grains/cell. Finally, a cell number titration and DIVESO titra- 
tion were completed. Similar statistics were achieved for all 
cell concentrations chosen (Fig. 6B), and DIMSO was found 
to have no effect at < 0.3% (Fig. 6C). Therefore, the screen- 
ing assay protocol was set to 4,000 cells per well with a final 
DIVISO concentration of 0.3%. 
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Fig. (5). Dose response to wortmannin and HC in Clone 10 cells. Each point represents Mean + SD for 6 replicates. HC inhibits the lysosomal 
degradation of the autophago somes so that more accumulate for imaging with increasing doses of HC. Wortmannin IC50 was similar at 5 and 
10 jiM HC (1.4 nM). A dose of 33 nM for wortmannin and 10 jiM for HC was selected to produce a good signal: background and best Z'. 





Fig. (6). Assay development. (A): Time of starvation of Clone 10 cells in the presence of wortmannin. Points show Mean + SD for 4 repli- 
cates. The signal was stable by 2 hr starvation and was not improved at 3 hr (Table). (B): Titration of number of Clone 10 cells/well, in the 
presence of wortmannin under starvation conditions for 2 hr. Bars show Mean + SD for 4 replicates. The response to wortmannin was similar 
at all cell densities and the Z' was stable at 0.80 + 0.04 across the range. A cell concentration of 4,000/well was selected. (C): DMSO titration 
for Clone 10 cells + 33 nM wortmannin under starvation conditions for 2 hr. Points are Mean + SD for 4 replicates. Table: for each DMSO %, 
the Z' was calculated from the 33 nM v^. 0 nM wortmannin value; Signal Window was the ratio of the same values. No effect of DMSO was 
seen < 0.3% DMSO. (D): Effect of binning images of Clone 10 cells, 4000 cells/well in the presence of wortmannin under starvation condi- 
tions for 2 hr. The same cell images were used with and without applying 2x2 binning. Points are Mean + SD for 4 replicates. No effect of 
binning on the Z' (0.82 v^. 0.79) or wortmann in IC50 values was observed. 
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Fig. (7). Assay Qualification. Left: Correlation plot in Spotfire from % inhibition of autophagy on Day 1 vs. Day 2 of the Validation Screen. 
From 8,000 compounds (25 plates) screened and setting a threshold of >30% inhibition of the positive control, the hit rate on Day 1 was 3.0% 
vs. 3.2% on Day 2. Confirmation of positive compounds between Day 1 and 2 averaged 69%). Right: Z' for control wells of individual plates 
(n=16 for each control, every plate) run for assay validation on Day 1 (pink) and 2 (blue). All plates were >0.5, with the Mean Z' being 0.66. 
These Z' were predictive of the HTS assay itself, indicating a robust assay. Where Z' <0.5, the plate was repeated. 



Image binning combines pixels together and sums their 
intensity. It increases scan and image analysis speeds and 
reduces the size of data files, and therefore is very beneficial 
for HTS. However, binning decreases image resolution. To 
determine if image binning was feasible for this assay, cells 
were treated as already described and imaged on the IN Cell 
3000. The images were then analyzed using either no bin- 
ning or 2x2 binning to determine % cells with > 6 grains per 
cell. Binning had no effect on Z', or wortmannin IC50 values 
(Fig. 6D) and so 2x2 image binning was adopted for primary 
screening. Under these conditions each 384-well plate took 
30 min to acquire and analyze. 

Primary Screening 

Transitioning from assay development to primary screen- 
ing involved an assay qualification step, whereby a "Valida- 
tion Library" of about 8,000 pre-selected compounds was 
run on two separate days using reagents and equipment con- 
figured identically with the planned screening process. This 
step allowed us to check the reproducibility of compound 
activity in the assay, get a projected expected hit rate, and 
also assess the plate to plate variability in terms of Z', before 
committing to an entire screen. A reasonable correlation was 
observed between % activities of compounds on the 2 days 
(Fig. 7A). From our assay development we were expecting a 
Z' for controls well above our minimum acceptable level of 
>0.4 for a cell-based screen, and this was achieved for all 
plates, with an average Z' of 0.66 (Fig. 7B). The assay was 
accepted into HTS and went on to test around 240,000 more 
compounds in the primary screen. The screening assay was 
run with the aid of a small robotics assembly consisting of a 



Thermo-Fisher CatX robot and plate "hotel" feeding a 
Combi and the plates were read off-line on the IN Cell 3000, 
usually overnight. Statistics for the primary screen followed 
those of the validation screen closely, with the excellent 
plate Z' being maintained and very few repeat assays re- 
quired. With a cut-off of >30% inhibition of grain formation, 
a hit rate of 2.8% was obtained, about 10,500 compounds 
(Fig. 8). 

Reanalysis of Archived Images 

As with many large compound libraries, a certain 
number of fluorescent compounds interfering with the GFP 
or Draq5 readout through actually binding to cells are pre- 
sent in the collection. Also, a number of compounds will be 
significantly cytotoxic even in a 2-hr assay. Example images 
of such compounds are shown in Fig. (9), A-F. Quantitative 
cell image analysis allowed us to set additional criteria to 
identify these categories. Therefore, following the comple- 
tion of primary screening, fluorescent, quenching and 
acutely toxic compounds were picked out by applying a sec- 
ondary data analysis. Highly fluorescent compounds staining 
cells are easily identified since they produce an increase in 
Icyt (average GFP intensity in cytoplasmic region), and this 
is often accompanied by an increase in Dmrk (average nu- 
clear diameter) because fluorescent compounds can interfere 
with nuclear mask determination. GFP quenchers produce a 
decrease in Icyt values. Toxic compounds can be identified 
using the features Npass, Icyt and Dmrk. Acutely toxic com- 
pounds reduce the number of cells analyzed per well 
(Npass), and compounds that do not reduce cell number but 
produce undesirable morphological changes such as nuclear 
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Fig. (8). Primary screen - Spotfire plots of results. Left : scatter plot (each point is a separate compound) with normahzed results for primary 
screening (% inhibition of autophagy vs. compound tested). Horizontal lines represent Mean with ranges + 2 x SD (Mean = 1.8 + 15.9). Com- 
pounds were selected positive at > 30% inhibition and are marked in red, rejected compounds in blue. Right : values binned as a histogram of 
activity distribution of all compounds marked as for the Scatter plot; inset is an expanded view of the distribution of selected positives. 





Fig. (9). Compound quahfication: example images (20x) of cells after treatment with compounds excluded through the secondary data analy- 
sis; Green - EGFP channel fluorescence, Red - nuclear channel fluorescence. (A): Positive control (0.3% DMSO); (B): Negative Control (33 
nM wortmannin); (C): Fluorescent compound; (D): EGFP quenching compound; (E): very toxic compound; (F): somewhat toxic compound. 
Data associated are shown in Table 2. (G): Re-analysis of primary screen image data to remove fluorescent compounds. Compounds were 
removed that increased the normalized Dmrk (average nuclear diameter) or Icyt (average intensity in cytoplasmic region) values by more than 
3 SDs from mean values (Mean = 0.06+6.05 - horizontal lines). 



shrinkage and cell rounding, will decrease Icyt and Dmrk 
values significantly as compared to controls. Such false posi- 
tive or undesirable compounds were filtered by eliminating 
those that reduce raw Npass values in 4 fields of view to 
below 200 cells (this represents about an 80% reduction in 



cell number) or increase, or decrease, normalized Dmrk and 
Icyt values by more than 3 standard deviations from mean 
control values (Fig. 9G). Example data from such com- 
pounds are shown in Table 2. 
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Table 2. Example Data for Filtering Compound Interference Using 4 Features 



Compound 


% Pos (Norm) 


Icyt (Norm) 


Dmrk (Norm) 


Npass (Raw) 


iNcgaiivc \^oiiiroi (yj.j /o uiviow ) 


-u.z 


1 f> 


n 7 

u. / 


Jy 1 


Positive Control (Wortmannin 33 nNI) 


101.3 


-5.4 


1.1 


577 


Fluorescent compound 


94.6 


56.1 


42.9 


631 


GFP quenching compound 


43.5 


-73.8 


-4.8 


627 


Very toxic compound 


94.5 


-48.3 


-26.2 


153 


Somewhat toxic compound 


107.2 


-46.3 


-17.0 


551 



'Norm' values are normalized to the mean of all plate wells; N Pass (Raw) is cell number in 4 fields of view 
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% inhibition of cell viability 



Fig. (10). Confirmation of the primary screen compounds (Y axis) vs. viability counterscreen (X axis). Compounds were run at 10 jiM in 
triplicate. Data are normalized as the % of controls on each test plate; each point is the Mean of 3 replicates. Unconfirmed screening hits are 
marked in blue; confirmed positives are marked in red. Compounds with > 50 % inhibition of cell viability in the counterscreen are in the 
green box. Of the confirmed positives, 122 showed > 50% inhibition of cell viabihty (red points in green box) and were eliminated from 
further consideration. 



Confirmation Assay 

After these filters had been applied, about 30% of the 
compounds were eliminated and a selection of 7,500 com- 
pounds were re-tested in triplicate at 10 |iM in a confirma- 
tion of the primary assay, with the same > 30% inhibition 
cut-off (Fig. 10). From this we confirmed about 2,000 (27%)) 
of the compounds. Interestingly, this low confirmation rate 
was directly brought about by our original rather "relaxed" 
hit criterion of >30%) inhibition: only 13%) of compounds 
with between 30 and 50% inhibition in the primary assay 
were confirmed, whereas 46%) of compounds with > 50%) 
inhibition confirmed (data not shown). 

Counterscreen 

The primary screen selected compounds that inhibited 
autophagy when the cells were placed in starvation medium. 



which would ultimately lead to cell death although not in the 
2 hr time-frame of the assay. Acutely toxic compounds that 
eliminated cells were removed in the secondary data analy- 
sis, but we also wished to eliminate compounds that were 
toxic to cells not undergoing autophagy and over a longer 
time period. Compounds successfully reaching the confirma- 
tion stage were therefore further tested for toxicity occurring 
in full growth medium, again at 10 |iM in triplicate, using 
the WT HeLa cell line growing in full growth medium for 24 
hr and a standard cell viability assay. Those compounds 
showing >50% toxicity (a further 122) were eliminated by 
this assay (Fig. 10). 

Chemical Triage 

At this stage, the remaining compounds were exam- 
ined for structural attributes and their previously documented 
activities in other screens. Those with undesirable properties 
(frequent hitters, non-specific inhibitors, undesirable chemi- 
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Compound ID 

Fig. (11). Distribution of IC50 values of -1,400 compounds tested in a dose response ranging from 2 nM to 30 jiM. About 1200 compounds 
gave IC50 values of < 30 jiM and of these, 280 compounds gave sub-jiM values (green); 79 compounds were classed as inactive (IC5o>30 
jiM) (orange). 



cal structure, etc.) were identified and around 350 were re- 
moved from consideration. 

Dose Response Determination 

Confirmed compounds (-1,400) were placed in a 10- 
point dose response (in triplicate) from 30 jiM to 1.5 nM. 
Those compounds that were demonstrably dose responsive 
and also achieved a maximum activity of >30% inhibition 
were considered active, and about 1,200 compounds fit these 
criteria. A range of activities were observed (Fig. 11) from 
nM to |iM. Nearly 300 compounds showed nM activity. 

Secondary Assay 

Finally, we added a longer-term cell death assay into 
the screening tree in a second cell type, the tumor cell line 
HI 299, to identify compounds that selectively brought about 
cell death under starvation conditions but not under normal 
growth conditions by 48 hr. From this about 400 compounds 
that were not toxic for either cell type in full nutrient condi- 
tions and only caused loss of viability under starvation con- 
ditions were identified (Fig. 12) and classified as "Active 
Specific". 

DISCUSSION 

Our HTS/HCA screening technology was effectively de- 
ployed to perform a good-sized screen of about 0.25 million 
compounds, and in hindsight, we could have done many 
more. The assay was simple once developed and both robust 
and relatively inexpensive. Compounds inhibiting EGFP- 



LC3 granule formation by >30% at 2 hr in starvation me- 
dium were designated positive (-7,000 out of 246,000 
tested). At the same time, parameters of cell health (number, 
nuclear size) could be collected to identify particularly toxic 
compounds; fluorescent compounds could also be identified, 
as well as GFP signal quenchers. An advantage of image- 
based assays is that such filters can be applied early and in a 
single assay, making these assays much "cleaner" for com- 
pound qualification. Compounds that were cytotoxic to WT 
HeLa cells given a longer exposure (24 hr) to compound, 
and when in full growth medium and not nutrient deprived, 
were also identified and removed from consideration. 

Applying these filters, moving from the primary screen- 
ing hit rate of 2.8% and around 10,500 compounds, after 
dose response only ~ 1,200 compounds (0.5%) remained, 
with IC50 values ranging from low nM to |iM in the EGFP- 
LC3 primary assay. For further compound characterization, 
both bio- and chemi- informatics were extensively employed 
to assign possible functions or cellular pathways to the com- 
pounds, based on both internal archived data for each com- 
pound or compound series, and on external literature. As 
examples (results not shown), we found compound series 
that inhibit PI3KC3/Vps34 using recombinant Vps34 protein 
in a TR-FRET format, and another family was known active 
against ULK from results with the Millipore Kinase Profil- 
er™ assay panel. We also found series for which the targets 
were unidentified. In order to fully qualify the compound 
activity in relation to the autophagy mechanism and also that 
activity was not confined to HeLa cells, we added a further 
assay pair into the screening tree in a second cell type, the 
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Fig. (12). Left : Dose response IC50 values (M) of -1,200 compounds which were designated Active in the Primary and Counter screen assays 
in EGFP-LC3-HeLa cells plotted vs. their cytotoxicity activity in the Secondary starvation assay in HI 299 cells. So that they appear on the 
plot, compounds are set to 100 jiM activity if inactive (red, n = -450), and those with IC50 values > 30 jiM are set to 50 jiM (n= -50). For the 
two secondary assays in HI 299 cells, compounds designated Active Specific (green) were compounds which gave IC50 values in the starva- 
tion medium which were 3x or more lower than in the parallel assay in fiill growth medium (-400 compounds). Active Non-specific com- 
pounds (yellow) have similar IC50 values in both formats. Right : Examples of compound dose responses obtained in the secondary assay 
(pink: starved; blue: full growth medium). (A): Active Non-specific, (B) and (C): Active Specific. 



tumor cell line HI 299, in full dose response mode, to select 
for compounds that only caused cell death under starvation 
conditions and not with full growth medium, by 48 hr. From 
this, about 400 compounds were identified as both active and 
specific for inhibiting starvation-induced autophagy. The 
final screening "tree" and outcome of compound triage is 
illustrated in Fig. (13). 

Different groups have already described that the two 
kinases Vps34 and ULK are key regulators of the autophagy 
machinery [1, 7-9]. The fact that we identified compounds 
active against these 2 potential drug targets for autophagy in 
our screen validates the design and the well-defined and 
planned path of our HTS on autophagy. Chemical matter 
needs to be further explored, including selectivity assess- 
ments and ADMET, in order to support discovery and move 
the project forward. Many autophagy inhibitors used in the 
autophagy field are in fact promiscuous compounds, such as 
3-methyladenine (3 -MA), or compounds with limited po- 
tency, such as hydroxychloroquine. Though there are already 
clinical trials ongoing using hydroxychloroquine for cancer, 
there remains an urgent need to identify potent and selective 
compounds for therapeutic applications. 

There has been a great deal of discussion in the field 
about the advantages and difficulties of "phenotypic screen- 
ing" although in reality, cellular assays have been routine in 
HTS for many years and any cellular assay (e.g. Ca^-rise 
readouts, cAMP readouts ...) is phenotypic until the in- 
tended molecular target is confirmed. While phenotypic 



screening has the apparent disadvantage of a lack of a known 
molecular target at the outset, the converse aspect is that 
setting up such assays does not require the prior selection of 
a molecular target which might turn out to be the wrong one. 
In addition, since these are cellular assays, compounds that 
cannot access intra-cellular targets or function in the pres- 
ence of physiological protein levels are weeded out early; 
and as screening libraries grow bigger and the quality better, 
the frequent though dubious argument that bad characteris- 
tics can be overcome later should become less employed. It 
should also be acknowledged that picking even an obvious 
molecular target has not historically created a straight path- 
way to drug discovery, and many approved molecules later 
reveal side effects or a polypharmacology that belies any 
specificity for the original target [10, 11]. An analysis of 75 
first- in-class small molecules approved by the FDA in the 10 
years before 2009, when the target-based HTS approach was 
at its zenith, showed that 28 came from phenotypic screening 
approaches where no target was defined beforehand, and 17 
from target-based approaches [12]. Deploying the huge 
chemical library resources of large pharmaceutical compa- 
nies in a phenotypic assay means that they are truly being 
assessed against the complex biological machinery that con- 
tributes to a cellular response in ways that might produce 
something entirely novel as well as biologically active. 

As we have ramped up our expertise in image-based 
screening and put the necessary data reduction methods in 
place, the screening of many more compounds has progres- 
sively become less of a technical challenge, with the next 
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Fig. (13). Screening tree of HTS process. 

major hurdle being the substitution into screening, for non- 
oncology targets at least, of "more relevant" cell types (e.g. 
primary cells; iPSC-derived cells). However, a further and 
harder to quantitate disadvantage that exists for a project 
completing a large-scale phenotypic screen is that it presents 
a generation of medicinal chemists that were trained largely 
on the target-based approach with a daunting amount of 
chemical data and a more obscure biology. After chemi- and 
bio-informatics triage of a compound set, it is often hard for 
project teams to distract attention from compounds that have 
a historic track record with a molecular target in favor of 
perhaps more novel pathways or as yet unknown targets in- 
teracting with other positive compound series or "single- 
tons". But enabling the required processes has been widely 
recognized as necessary [13, 14] and tools are increasingly 
available [15, 16]. Given the striking success of phenotypic 
approaches to deliver first- in- class therapies [12], devising 
and implementing higher throughput approaches to post- 
screening, biomedical/translational qualification of com- 
pounds should be the next frontier of the continuing chal- 
lenge to deliver new medicines for unmet medical needs. 
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ABBREVIATIONS 



HTS = High Throughput Screening 

HCA = High Content Analysis 

LC3 = microtubule- associated protein 1 light-chain 

3 (MAP1LC3) 

Atg = Autophagy related gene 

GFP = Green Fluorescent Protein 
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BBSS 


= Earle's Balanced Salt Solution 


DPBS 




HC 


= Hydroxychloroquine sulfate 


DMSO 


= Dimethyl sulfoxide 


ADMET 


= Absorption Distribution Metabolism Excre- 




tion Toxicity 
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